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Abstract 

We suggest that the mass of four compact stars SAX J1808.4— 3658, KS 
1731-260, SAX J1750.8-2900 and IGR J17191-2821 can be determined 
from the difference in the observed kiloHertz quasi periodic oscillations (kHz 
QPO-s) of these stars. The stellar radius is very close to the marginally stable 
orbit R ms as predicted by Einstein's general relativity. It may be noted that 
the first of these stars was suggested to be a strange star more than a decade 
back by Li et al. (1999) from the unique millisecond X-ray pulsations with 
an accurate determination of its rotation period. It showed kHz QPO-s eight 
years back and so far it is the only set that has been observed. This is the 
first time we give an estimate of the mass of the star and of three other 
compact stars in Low-Mass X-ray Binaries using their observed kHz QPO-s. 
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1. Introduction 

Strong variability on millisecond time scales, the kilohertz quasi-periodic 
oscillations (kHz-QPOs) observed in the low mass X-ray binaries (LMXBs) 
was first reported by van der Klis et al. (1996), with observations from the 
Rossi X-ray Timing Explorer (RXTE). Ever since, it led to many theoretical 
modelling and prediction of its origin, and to determine the correct equation 
of state (EOS) of the compact objects (see van der Klis, 2006 for a recent re- 
view). Although in its earlier stage of detection it was found in the Atoll and 
Z class of sources of the LMXBs, later it was also observed in the millisecond 
X-ray pulsars and a few uncategorised sources. 

Typically these characteristic frequencies have two peaks, the lower iyiow) 
one ranges from a few hundreds of hertz and the upper (v up ) one often goes 
beyond a kilohertz. Although some QPOs with both the frequencies falling 
much below the kilohertz range have been observed in some black hole can- 
didates, for the present study we shall refrain from considering them. 

Here we start from a recent argument if the kHz-QPO-s are related in 
any sort of way to the spin frequency of the star. Mendez and Belloni (2007) 
re-examined all reliable available data of sources for which there exist mea- 
surements of two simultaneous kHz QPO-s and spin frequencies (v s ) and 
advance the possibility that the difference between the QPO-s, Av (the dif- 
ference between the upper and the lower frequencies) are not related to each 
other. We have used their data for these ten sources and we find a correlation 
between the two simultaneous kHz QPO-s raised to the power 2/3 (Fig [I]). 
Since QPO-s are represented by peaks in the power spectrum, they are re- 
lated to the relatively stable orbits of accreting matter from which particles 
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Figure 1: Plot of v up vs. v^ w for different LMXBs. Different symbols and colours have 
been used for different sources. 



fall onto the compact object - in our case assumed to be a strange star - 
and these orbits are the maximally stable orbit given by an average radius 
of Rms for h> up and the outermost orbit Ro for V[ ow . This assumption leads to 
possible determination of the mass and radius of the compact strange star. 

Turning to history, the most popular model to explain the difference of 
kHz QPOs was the beat frequency model (Strohmayer et ai, 1996) which 
assumed the upper QPO (u up ) as the Keplerian orbital frequency Uk(r) of 
the innermost orbit in the accretion disk around the star, the separation 
between the two peaks i.e. Au peak = v up — u [ow as u spin and the lower QPO 
iyimu) as the beat frequency of v\.{r) and v svin . This model had suggested 
that for a particular star Au = v sp i n will remain constant. However, it has 
been observed that in many sources the Au is not related to z/ sp j ra .Also Au 
changes with time. The QPO resonance model of Kluzniak and Abramowicz 
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(2001) is a beautiful model but comes with its own set of problems. It 
states that the QPO could be explained by non linear resonant motions of 
accreting fluid in strong gravity. The two QPO frequencies depend here on 
the amplitude of the oscillatory motion and are approximately in the ratio 
3:2. However, the difficulties mainly arise from the fact that QPO frequencies 
are not consistently equal to this mentioned ratio, and the twin peaks shift 
their positions (Rebusco, 2008). Although all of the present QPO models 
come with some conceptual difficulty or the other, for our present work, we 
shall adhere to the QPO model of Titarchuk and Osherovich (1999) which 
allowed us in an earlier paper to determine the mass and limits to the radius 
of the accreting star (Li et al. 1999). 

Most compact stars showing a pair of QPO-s in the X-ray power spectrum, 
have variable 5v. The first star to show accurately constant X-ray pulse 
period was the X-ray pulsar SAX J1808.4— 3658. This source was found 
to show kHz QPO-s with the a frequency difference of 195 kHz. This is 
roughly half the spin frequency of the star u spin which is measured very 
accurately to be 401 Hz. So a modified beat frequency model was tried with 
Vpeak ~ 0.5z/ sp j„. But for many other stars the difference of kHz QPO-s has 
no apparent correlation with the spin frequency of the star like that has been 
suggested by the beat frequency model. 

There are now at least five other X-ray pulsars whose spin is steady and 
accurately measured and only for one of these, XTE J1807— 294 kHz QPO-s 
have been seen. We find, from applying the Titarchuck and Osherovich(1999) 
model, that the masses (M) and radii (R) of these two stars are very small 
and may fit those of strange stars. 
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2. Details of Calculation and Results 

2.1. Derivation of the mass and a limit to the radius of the compact star, 
using QPO frequencies 

We assume the possibility that the lower QPO frequency is due to parti- 
cles thickly clustered at the inner edge of the accretion disc (at a radial dis- 
tance R ) and the upper QPO frequency is due to particles thickly clustered 
at the marginally stable orbit (at a radial distance R ms )- The expression for 
the marginally stable radius (Bardeen, Press and Teukolsky, 1972) is given 
as: 

R ms = R g {3 + Z 2 - [(3 - Z0(3 + Z x + 2Z2)} 1 ' 2 } (1) 

where 

Z x = 1 + [1 - (a/i? 9 ) 2 ] 1/3 [(l + a/R,) 1 / 3 + (1 - a/R g )^] (2) 

Z 2 = [^{a/R g ) 2 + Z 2 ]^ (3) 

where R s = 2GM/c 2 , R g = GM/c 2 and a = 2it v &V i n l /Mc, I is the moment 
of inertia, c is the velocity of light and G is the gravitational constant. A 
lower limit for R according to the Titarchuck and Osherovich model is as 
follows (Li et ai, 1999b) 

/ M \ 1/3 

«o = 3fl,x(-j (4) 

and hence it is clear that i? depends only on M and whereas R ms depends 
on both M and v svin . So for a particular star (fixed M and v svin ) R ms should 
be constant. From the Fig((T]) above we see that vlp 3 /uf^ ~ 1.2 — 1.3 for 
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most of the stars and one can calculate the most probable outermost orbit 
r and consider its relation to the Ro given above. We tabulate vi ow , v up and 
Uspin for 10 sources in Tabled] for which Fig([T]) shows that Vup frf^ is nearly 
constant. 



Table 1: kHz QPO frequencies and spin frequencies for 
different stars 



Sources 


^low 


L'up 


^ 'spin H Z 


X-Ray pulsars 


SAX J1808.4-3658 


499 


694 


401 


XTE J1807-294 


127.1 


352.8 


191.0 




184 


374 






195 


393.4 






216.5 


449.6 






246.3 


466.3 






253.9 


490.6 






359.9 


563.6 




other stars 


KS 1731-260 


898.3 


1158.6 


524 


IGR J17191— 2821 


900 


1227 


294 


SAX J1750.8-2900 


936 


1253 


600.75 


4U 1728-34 


576.6 


925.0 


363 




614.9 


967.6 






647.1 


995.2 






702.3 


1053 
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749.7 


1094.2 






793.0 


1148.9 






853.8 


1183.2 






894.2 


1172.9 






308 


582 




4U 1608-52 


472.8 


798.3 


619.0 




555.6 


872.5 






561.3 


874.5 






587.3 


893.0 






623.2 


922.7 






660.0 


956.5 






682.5 


993.7 






725.0 


1040.9 






754.0 


1044.9 






772.3 


1024.9 






787.7 


1047.0 






822.7 


1061.5 






866.9 


1092.2 




4U 1636-53 


644 


971 


582 




688 


1013 






723 


1053 






769 


1058 






841 


1102 






856 


1114 






871 


1122 







886 


1126 






901 


1150 






913 


1156 




4U 1702-43 


657.1 


1000.1 


330 




707.7 


1037.7 






770.0 


1084.8 




4U 1915-05 


706.9 


1055.3 


272 




656.9 


906.4 






558.3 


911.2 






583.2 


923.7 





From equations (1) to (4) we obtain the difference 

v up — vi ow = [(R /R ms ) 3 ^ 2 — l]viow, (5) 

so that we can seek the stellar mass, the corresponding stellar radius and the 
moment of inertia to fit this difference. If the difference remains constant 
as in the case of SAX J1808.4— 3658 (at 195 Hz) we can predict the mass 
to be 0.904 M . For this mass and spin the R ms = 7.409 km. For strange 
stars, the corresponding radius is known from our EOS (see section 2.2) to 
be 6.89 km which is nicely inside R ms . The derived value of r = 8.607 
km and the approximate relation R = 3 x 2.9532M 1 / 3 gives a result 
R = 8.566 km which is very close to r such that our model leads to a 
consistent result. A recent estimate of the mass of SAX J1808.4— 3658 by 
Elebert et al, 2009, from its optical spectroscopy and photometry showed 
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that its mass is 0.9 ± 0.4 M Q . This is consistent with the mass derived 
according to the strange-star model. 

2.2. The mass-radius relation of strange stars 

For use in this paper, we have fitted the mass M/M & of the strange star 
model (Dey et al, 1998) along with the corresponding moment of inertia /, 
as a polynomial of the radius R, for the range where the radius increases 
monotonically with mass as follows: 

M/M Q = alR bl + a2R b2 (6) 

/ = glR hl + g2R h2 (7) 

The parameters are given in Table [2j This is for the EOS A which gives 
the right masses and magnetic moments of baryons in a finite system calcu- 
lation with our chosen interquark force as given in Bagchi et al. (2006). We 
would like to stress that this SS model has the Richardson potential modi- 
fied with two appropriate but different scales for confinement and asymptotic 
freedom for the interquark force. The parameters of this force are all fixed 
from elaborate baryon mass and magnetic moment calculations and satisfies 
the stringent binding energy requirement that 

E/A(uds) < E/A(Fe 56 ) < E/A{ud). (8) 

2.3. Application to other Low-Mass X-ray Binaries 

For KS 1731-260, SAX J1750.8-2900 and IGR J17191-2821, the masses 
and radii are 0.878, 0.895, 0.843 in M & and 6.84, 6.88, 6.762, R ms = 7.00, 
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Tabic 2: Parameters for fitting mass and moment of inertia for EOS A 



Parameters 


M/M 


Parameters 


I 


al 


4.6 x 1(T 2 


gl 


8.9 x 10" 4 


a2 


3.0 x 10- 6 


g2 


1.7 x 10~ 7 


bl 


3.0 


hi 


5.0 


b2 


7.2 


h2 


9.0 



7.031, 7.051 in km respectively. One can see the closeness of R ms to the stellar 
radius and this was what had made the estimate for 4U 1728—34 reliable in 
the paper by Li et al, (1999b). 

For the other millisecond pulsar XTE J1807— 294 if we may estimate the 
mass and the radius to be M = 0.782 M and R = 6.62 km with the 
R ms very close ~ 6.672 km and Rq to be 8.162 km - then for the seven sets 
of QPO-s given by Mendez and Belloni (2007), we find the r to be 13.177, 
10.705, 10.652, 10.859, 10.210, 10.350 and 8.997 km. Since R is just a limit 
- our results could mean that the outer orbit varies for this and the other 
stars for example, due to variations in the accretion rate. 

3. Discussions 

For the last few years, various people have searched for a relationship 
between QPO frequencies and spin frequency of the star. As an example, 
Belloni, Mendez and Homan (2005, 2007) have shown that v\ ow and v up are 
strongly correlated along a curve passing through the origin. Again, Yin et 
al. (2007) showed that, for six Atoll sources, the relation between Au and 
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v sp in can be fitted by a linear relation 

(Au) = -(0.19 ± 0m)u spin + (389.40 ± 21.67) Hz (9) 

However, according to these authors, this anti-correlation is just a conjecture 
since it is based on data of only six sources. As we have stated before, 
examining all available data of sources for which there exist measurements 
of two simultaneous kHz QPO-s and u s , Mendez and Belloni (2007) have 
advanced the possibility that Au and v s are not related to each other. On the 
other hand, we find a correlation between Uup and i/f^ which is expected from 
the Titarchuck and Osherovich model (our preferred model) in agreement 
with the findings of Belloni, Mendez and Homan (2007). Though presently we 
have data for ten sources only, we hope with future advanced X-ray satellites, 
like the ASTROSAT satellite to be launched by India in 2011, we shall be 
able to get more data on kHz QPO-s which will help us to re-examine findings 
of the present paper. 
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